Integrated optimization of scheduling for unmanned follow-me cars on airport surface

To promote the application of automated vehicles in large airports, in this study, we present an integrated optimization method for scheduling Unmanned follow-me cars. The scheduling process is divided into three phases: Dispatch, Guidance, and Recycle. For the Dispatch phase, we establish a vehicle assignment model, to allocate the vehicle resource equitably. For the Guidance phase, we offer an quantitative way, to measure the spacing between Unmanned follow-me car and aircraft. To optimize the efficiency of airport operation in the three phases and ensure safety, the collaborative planning model, and the conflict prediction model are established. An improved grey wolf optimization algorithm is adopted to enhance the convergence speed and generalization performance. A case study at Ezhou Huahu Airport in China demonstrates the effectiveness of the methods. The results show that the model of collaborative planning can make the balance of path selection, Unmanned follow-me car’s working time, and departure sequence. The convergence speed of the improved algorithm has been increased by 18.75%. The inequity index of vehicle assignment is only 0.015731, and the spatiotemporal distribution of conflicts is influenced by the airport’s surface layout.


Integrated optimization of scheduling for unmanned follow-me cars on airport surface Dezhou Yuan , Xinping Zhu * , Yajun Zou & Qing Zhao
To promote the application of automated vehicles in large airports, in this study, we present an integrated optimization method for scheduling Unmanned follow-me cars.The scheduling process is divided into three phases: Dispatch, Guidance, and Recycle.For the Dispatch phase, we establish a vehicle assignment model, to allocate the vehicle resource equitably.For the Guidance phase, we offer an quantitative way, to measure the spacing between Unmanned follow-me car and aircraft.To optimize the efficiency of airport operation in the three phases and ensure safety, the collaborative planning model, and the conflict prediction model are established.An improved grey wolf optimization algorithm is adopted to enhance the convergence speed and generalization performance.A case study at Ezhou Huahu Airport in China demonstrates the effectiveness of the methods.The results show that the model of collaborative planning can make the balance of path selection, Unmanned follow-me car's working time, and departure sequence.The convergence speed of the improved algorithm has been increased by 18.75%.The inequity index of vehicle assignment is only 0.015731, and the spatiotemporal distribution of conflicts is influenced by the airport's surface layout.
With the recovery of civil aviation after the epidemic of COVID-19, the introduction of unmanned vehicles at airports is considered as an innovative measure to improve airport operational efficiency and reduce costs.These vehicles, equipped with advanced sensors, controllers, enable autonomous driving functions with complex environment perception, intelligent decision-making, and collaborative intelligence.They have become the new generation of airport equipment for operations and support 1 .However, airport operation is a critical issue about safety, and it requires full consideration.Now that a single unmanned vehicle, has the capability for automatic control 2 , the overall risk control of the aviation transportation system in the hermetic environment of airport surface has become the focus.How to conduct scheduling of the fleet scientifically is a major problem for airports.
To balance the security and efficiency of the fleet, scholars have conducted research on the scheduling problem of unmanned vehicles in the airport scene.Basically, the scheduling process involves path planning and vehicle assignment 3 .For path planning, an optimized combination model is introduced.Because this model needs to discuss the situation in different velocities of the vehicles, it was solved with the heuristic, such as hill climbing particle swarm optimization (HC-PSO) algorithm 4 .For the other, the utilization ratio is considered as a kind of index, to evaluate the equity of vehicle assignment.A dynamic programming model with various charging strategies meets the needs of allocation for vehicle resource 5 .Generally, the scheduling problem is a classic vehicle routing problem with time windows (VRPTW).Different models and algorithms have been explored to address this challenge, including linear programming (LP), combinatorial optimization, graph theory, and network analysis 6 .Common algorithms include exhaustive methods, C-W savings algorithms, column generation algorithms, heuristic search, genetic algorithms (GA), and particle swarm optimization (PSO) algorithms.However, research on unmanned vehicle scheduling on airport surface mostly focused on electric tractors, which are subject to traction speed limitations 7 and do not meet the requirements of large-scale airports.There is a lack of research on the scheduling process of unmanned follow-me cars (UFMCs).
Regarding VRPTW on airport surface, integer linear programming (ILP) algorithms is utilized to determine the optimal allocation of ground support equipment (GSE) based on cost and time constraints 8 .The sorting of flight support operations for airport service vehicles is improved by GA, to optimize the structure and quantity of service vehicles.The objective is to enhance the efficiency of service vehicle utilization 9 .An energy consumption model for GSE is developed by the adaptive large neighborhood search (ALNS) algorithm, and provides scheduling solutions for scenarios involving both fuel-powered and electric vehicles, addressing the challenge of minimizing energy consumption 10 .Dispatching issue for electric GSE in the apron is addressed by hybrid neural networks (HNN), contributing to the overall improvement in efficiency 11 .An intelligent dispatching model for support vehicles is solved by GA, with a focus on optimizing the allocation of GSE and efficiency 12 .A theoretical system based on ILP is developed to solve the assignment problem of airport operational units, providing an effective framework for optimizing the allocation of operational units within the airport 13 .Collaborative decisionmaking for parking allocation and shuttle scheduling is successfully implemented through column generation (CG) algorithm, aiming to improve the efficiency of parking allocation and shuttle scheduling processes 14 .An objective function that minimizes operating costs and flight delay costs and optimizes GSE scheduling is devised with the help of ILP, to reduce operational costs and minimize flight delays through GSE scheduling 15 .Heuristic search methods are employed in a real-time scheduling decision system to efficiently assign GSE 16 .
Regarding aircraft path planning on airport surfaces, an ILP mathematical model for conflict-free aircraft trajectories in various traffic density scenarios is established with a rolling time domain calculation strategy 17 .A mixed ILP (MILP) model for gate allocation is developed by combining operation modes of multiple runway airports, aiming to optimize gate allocation 18 .Heuristic methods are implemented to search for the optimal speed profile for ground taxiing process, providing valuable references for aerodrome control units 19 .Optimal taxiing speed profiles were generated with an exhaustive method based on required time of arrival (RTA) 20 .Surface routing planning problems are also studied through multi-agent simulation 21 .
In conclusion, there are various optimization techniques employed in scheduling UFMCs on airport surfaces, including ILP, GA, ALNS, HNN, and heuristic search methods.Additionally, mathematical models such as ILP and MILP are utilized for path planning on airport surface.Some professional concepts such as RTA, energy management and flow of transportation management are introduced, that makes mathematical models meet the demands of business processes.However, further exploration is needed to develop more efficient and accurate optimization techniques for the scheduling of UFMCs on airport surface.
While previous research has provided insights into the issue of scheduling on airport surfaces, these studies have some limitations.Firstly, they mainly consider vehicles for small to medium-scale operations, typically involving up to 50 flights, which may not meet the demands of larger-scale vehicle scheduling.Secondly, although the vehicle assignment models may generate satisfactory solutions, they often fail to take uncertainties and flight support requirements into account.Finally, the evaluation of potential movement conflicts is not thoroughly addressed in these studies.The conflicts are either treated as disturbances or some mitigation suggestions are provided.
On the other hand, in advanced surface movement guidance and control systems (A-SMGCS) Level 4 and above, the surface movement of aircraft relies on three functional modules: target surveillance, conflict alert, and routing.However, the conventional method of surface taxiing guidance with lights of centerline on the taxiway proves to be expensive in terms of maintenance.It is not suitable for large and busy airports.To address this issue, an alternative approach within a generalized A-SMGCS implementation framework is proposed, which utilizes UFMCs as a replacement for taxiing guidance.Figure 1 shows the basic working principle of this approach, where the UFMC scheduling system integrates with the existing functional modules.In this concept, the aircraft follows the UFMC during taxiing, requiring close attention from the pilot for guidance.The advantage of UFMCs over traditional manual follow-me cars lies in their ability to formulate efficient vehicle scheduling schemes, enabling prompt departures from the parking lot and seamless integration with guidance tasks.Furthermore, UFMCs overcome the drawback of slower towing speeds associated with electric tractors, ensuring the efficiency of Figure 1.Working principle of the UFMC scheduling system.airport surface operations.Therefore, this paper aims to study the integrated optimization of UFMC scheduling to address this guidance requirement.

Problem description
The integrated scheduling process for UFMCs to facilitate aircraft taxiing guidance is illustrated in Fig. 2.This process comprises vehicle trajectory planning and vehicle assignment.Vehicle trajectory planning involves generating optimal trajectories for the UFMCs across three phases: dispatch from the parking lot to the start of guidance (Dispatch), implementation of the guidance process (Guidance), and return to the parking lot after completing the guidance (Recycle).These trajectories are generated based on specific safety and efficiency objectives.Each vehicle's trajectory consists of a sets of position coordinates x, y along the route, corresponding velocities (v) , and timestamps (t) .Vehicle assignment refers to the selection of a specific UFMC from a pool of candidates based on certain optimization objectives to perform ground taxiing guidance for the corresponding flight.
As can be observed, the precise planning of UFMCs in terms of starting, stopping, and speed profiles enables a higher level of refinement compared to manually operated follow-me cars.This allows vehicles to promptly reach the starting point of the guidance task and proactively plan to avoid potential conflicts during their operations.For the convenience of discussion, the UFMC and the aircraft being guided by it are considered as a "Guidance Unit (GU)" in the surface movement scenario.
From Fig. 3, during the Guidance phase, there exists a following relationship between the aircraft and the UFMC.The spacing between them affects the spatial occupancy of the GU.It also influences the scheduling decisions.To address this, a definition and a quantification method for the distance are provided: Definition 1 (Guidance following spacing) Within the GU, the spacing L from the front of the aircraft, to the rear of the UFMC.
As shown in Fig. 4a, the measurement of L is similar to the concept of Following Spacing in transportation engineering 22 : (1)  where l A and l V are the lengths of the aircraft and the UFMC, respectively.d AV represents the distance between their geometric centers.
During the following process, the visual influence to the pilot must be considered in order to comprehensively evaluate the safe guidance following spacing for each type of aircraft.As depicted in Fig. 4b, there exists the following relationship between the visual spacing d visual in the cockpit, and the guidance following spacing L is calculated by the Pythagorean theorem.

Model formulation
The scheduling model for UFMCs consists of two modules: vehicle trajectory planning and vehicle assignment.Figure 5 shows the details about this model: • The vehicle trajectory planning module comprises the method for guidance trajectory deduction and a conflict prediction model.The former is responsible for rapidly deducting and generating a solution set of trajectories based on pre-planned vehicle routes and delivery times, with the input of guidance tasks.The latter is used to evaluate potential conflicts and determine the feasibility of the solution set.Finally, all feasible solution sets are integrated to generate an optimal solution set of trajectories.• The vehicle assignment module transforms the optimal into corresponding guidance tasks, employing the UFMC assignment model to achieve the assignment function.Ultimately, the optimal scheduling scheme, i.e., the timetable of each guidance task, is obtained.

Collaborative planning model for surface guidance trajectories
The set of flights to be guided denoted as G = g 1 , g 2 , . . ., g s , where s means the total number of guidance tasks, i.e., the total number of flights.The three phases of the guidance task are denoted as P = p 1 , p 2 , p 3 , corresponding to the Dispatch phase, Guidance phase, and Recycle phase.

Objective function
To meet safety and efficiency requirements, the objective function is designed from the following perspectives: • Minimizing the number of conflicts The number of conflicts be denoted as N conflict , and set the objective to minimize the number of conflicts: (2) min N Conflict  www.nature.com/scientificreports/ • Minimizing the working time of UFMCs Minimizing the working time of UFMCs means reasonable route planning, and reducing the ineffective waiting time of UFMCs at the gap between two phases as possible.Therefore, the objective is to minimize the working time of UFMCs: Here, i represents the i-th guidance task in the timetable, j represents the phase at which the UFMC is located, and s is the total number of flights in the timetable.The working time of UFMCs is the sum of the working times in each phase of the guidance tasks, and it is measured in seconds.Therefore, it could reach the magnitude of 100,000.Meanwhile, the number of conflicts is up to 1000 extremely.
To achieve the comprehensive optimization of the two objectives, the range of working time of UFMCs is normalized to [0,1] as follows: Here, the working time of UFMCs is transformed into T ′ , to remove the influence of dimension while ensuring the monotonicity of the indicator.
Therefore, the objective function for vehicle trajectory planning is: While setting θ 1 = 1 , the integral part of the objective function value Z maps the number of conflicts on the field, and the decimal part maps the total working time of UFMCs on airport surface.The value of θ 2 is determined by the scale of the airport surface, ensuring that ∀N conflict ∈ N , θ 2 N conflict is one order of magnitude larger than θ 1 T ′23 , achieving simultaneous optimization of the two components.

Constraints
• Constraints of routing for UFMCs For the three phases of vehicle guidance, the UFMC scheduling system has prestored one shortest route and (k * − 1) sub-optimal routes for each pair of origins and destinations.During the process of scheduling, the system selects one from the k * candidates as the route of UFMC, i.e., where R ij is an integer decision variable that represents the number of the route.The subscript i and j indicate that route selection is performed for each guidance task and each phase in the trajectory planning module.

• Constraints of time windows for UFMCs
The aerodrome control unit and the Airport Operation Center (AOC) adjust the delivery time windows of the UFMCs through collaborative decision-making throughout the flight process to achieve integrated optimization of scheduling of UFMCs and control of flights: where γ ij is the decision variable, corresponding to the delivery time adjustment for the i-th guidance task in phase j. δ j represents the upper limit of adjustment for each phase.Taking guiding flight i as an example, the delivery times for the three phases of the UFMC are shown in Table 1.ETA and ETD are obtained from the flight schedule.They represent the Estimated Time of Arrival and Estimated Time of Departure, respectively.LDR is the time from the arrival aircraft touching down, to reaching the starting point of the Guidance.LUP is the time from the departure aircraft waiting at the runway entrance to having clearance of lining up the runway, respectively.CRS represents the time taken for the UFMC to pass through the taxiway on the apron, while UT is the time for the UFMC to make a U-turn in front of an arrival aircraft.
• Constraints of domain To ensure the meaningfulness of each variable, it is necessary to define their domain, including the number of conflicts as a natural number, the route number of the UFMC as a positive integer, and the non-negativity of the UFMC working time and delivery time adjustment:

Trajectory deduction
The purpose of guidance trajectory deduction is to deduce the trajectory of an UFMC or GU, for each guidance task at different phases based on a predetermined route and the delivery time.For UFMCs, precise control to parameters can be achieved based on the setting of guidance speeds in different zones, vehicle performance parameters, as well as road conditions, enabling accurate trajectory deduction to support the prediction of potential conflict in trajectories.

Definition 2 Trajectory of UFMC
The coordinates the UFMC's geometric center passes through, along with corresponding timestamps and velocities.

Definition 3 Trajectory of GU
The coordinates the GU's geometric center passes through, along with corresponding timestamps and velocities.
Since the preplanning process has generated several candidate routes, in order to accurately match the UFMC with the guidance requirements, further algorithm design is required, to generate high-precision velocity profiles and achieve synchronization between the actual delivery time and the start of the timestamp.
The input to the algorithm of Speed Profile Generation for UFMC is a set of nodes formed by a single candidate route.Then, correct the speed limit and determine the safety throttle/braking distance for each road segment, by basic kinematic equations.Finally, generate the final velocity profile: The algorithm is showed as following: Table 1.Delivery time of UFMC guiding flight i in three phases.

Flight Type
Phase Regular delivery time (RDT) Actual delivery time (ADT) Algorithm 1. Speed Profile Generation for Unmanned Follow-me Car.
The process for the guidance trajectory deduction is as follows: • Input the current phase j, the number of parking stand where flight i is located, and the selected route R ij .Generate the set of nodes x 0 , y 0 , x 1 , y 1 , x 2 , y 2 , . . ., x t , y t corresponding to the guidance route; • Input the delivery time adjustment γ ij and determine the delivery time, denoted as t 0 , according to Table 1 www.nature.com/scientificreports/Where v node represents the speed limit at the next node, i.e., the current section's speed limit.E Apron , E Taxiway , and E Turn represent the apron (including service roads), the taxiways in the maneuvering area, and the turn- ing points, respectively; • Call Algorithm 1 to generate the velocity profile {(v 1 , t 1 ), (v 2 , t 2 ), . . ., (v t , t t )}: • Smooth the segments at the turning areas as shown in Fig. 6. ∀node ∈ E Turn , adjust the timing as follows: Where α is the turning angle, r is the turning radius, x is the advance distance for the turn, s is half of the arc length corresponding to the turning angle, and t ′ Turn is the adjusted timing for passing the turning point.• Output the guidance trajectory:T i, j = x n , y n , v n , t n |n = 0, 1, 2, . . ., t .
The guidance trajectory T i, j provides information about the working time of UFMC: Moreover, the guidance trajectory T reflects the situation of the UFMC on surface, and serves as the foundation for the conflict prediction model.

Conflict prediction
A conflict prediction model based on protected zones is established by using the coordinates of each position in the trajectory as the geometric centers.The properties of a protected zone are determined by trajectory's parameters of corresponding position, enabling the prediction of conflicts.Among commonly used protected zones, the elliptical demonstrates good prediction accuracy and intuitiveness.
The lengths of the semi-major axis "a" and the semi-minor axis "b" of the elliptical protected zone are calculated as following: Here, k sec urity represents the safety factor.l target and w target denote the length and width of the UFMC or GU, respectively.During the Guidance phase, the moving target refers to the GU.During the Dispatch or Recycle phase, the target refers to the UFMC: Here, w A represents the wingspan of the guided aircraft, and w V represents the width of the UFMC.
The purpose of setting the safety factor k sec urity is to minimize the possibilities of both false alarms and missed alerts in conflict prediction.Figure 7 shows the scenarios of different safety factors corresponding to the protected area.
The black protected zone represents a safety factor of 1, which cannot envelop the entire moving target most of the time, potentially leading to missed alerts.The red protected zone corresponds to a safety factor of 1.2, which (15) can envelop the entire moving target and minimize the probability of false alarms when potential conflicts are predicted.The Reference 24 has already discussed incidents of wingtip collisions on the apron, it is sufficient for their model to recognize and deal with such incidents.Therefore, our model eliminates the wingtip clearance ε proposed by Giuseppe S [7], and uses the length of the GU instead of the length of the aircraft in their tractor problem to ensure that there is sufficient space for our GU within the protected zone.Figure 8 visually shows a typical scenario of conflict between GUs.The elliptical protected zone could accurately represent the spatial occupation of GU, and maintain safe spacing both longitudinally and laterally.To reproduce the protected zone, and determine whether the ellipses overlap or not, We have formulated the elliptical equation as following: Subsequently, the Algorithm 2 is used to output the number of conflicts, denoted as N Conflict , thereby completing the conflict prediction.

Model for UFMC assignment
Compared to traditional manually-operated vehicles that rely on human experience for scheduling, operation, and maintenance, UFMCs overcome the limitations of manual driving through fully automated control.Building upon precise trajectory planning, further considerations for vehicle maintenance and optimal allocation of vehicle resources are necessary.In this regard, an UFMC assignment model shown in Fig. 9 is designed in the UFMC scheduling system, wherein the optimal solution set of trajectories is transformed into specific guidance tasks, to achieved the function of assignment through ILP.www.nature.com/scientificreports/Considering that the source assignment problem can attribute to 0-1 programming.According to the general method in management of civil aviation 13 , we designed a model for vehicle assignment as following:

Model assumptions
• The UFMCs at the end of the Recycle phase must return to the parking lot where they started, in order to perform the next guidance task.• A buffer time of η is allocated between two guidance tasks to ensure orderly traffic flow in and out of the parking lot, or transitioning to the corresponding parking lot for the next guidance task.• A vehicle work continuously until its battery level reaches a low state, then return to the parking lot for a full recharge.The charging time C and the safe endurance q are fixed values specific to each type of vehicle.Factors such as battery decay and seasonal variations are ignored, but the buffer time η is included in the continuous working time.
Let F = f 1 , f 2 , . . ., f k represent the set of k UFMCs in the airport.

Objective function
This model aims to achieve the most equitable vehicle assignment, so we adopt the Gini coefficient as the evaluation metric.In the field of economics, a Gini coefficient between 0.3 and 0.4 is considered fair 25 .However, www.nature.com/scientificreports/considering that the working time of UFMCs is a "negative asset", the Gini coefficient is referred to as the "assignment inequality index" in this context.Therefore, the objective function is set to minimize the assignment inequality index of the working time for UFMCs: Where P l = 1 k , and Y l represents the proportion of working time for the l-th vehicle to the total working time, and V l is the cumulative distribution of Y l . 3j=1 t ij denotes the accumulated time for the three phases of a single guidance task.x il is a binary decision variable indicating whether vehicle l guides flight i.

• Constraints of flight coverage
Each flight should be assigned to one UFMC ONLY:

• Constraints of consecutive guidance tasks
The time windows of guidance tasks assigned to the same vehicle should not overlap.Let o and p represent a pair of guidance tasks corresponding to flights where their time windows overlap each others, then the constraints can be formulated as follows: The condition that the time windows of guidance tasks o and p overlap can be expressed as follows: Here, a o and a p correspond to the actual delivery times (ADT) of guidance tasks for flights o and p, respectively.
Usually, guidance tasks are sorted by ADT.Fig. 10a shows the overlap situation of time windows between adjacent index of guidance tasks.
• Constraints of vehicle electricity (22)    www.nature.com/scientificreports/Constraints of vehicle electricity ensure that the quantity of electricity of the UFMCs meets both upper and lower limit requirements during operation and forms a closed loop for charging: Here, d il is the binary representing whether vehicle l requires charging after completing task i. q reserve is a backup of the quantity of electricity that ensures the UFMC can safely return to the parking lot after operating for a continuous duration of q.
• Constraints of charging For single vehicle, the guidance task should not have conflict with the charging process: The conflict condition between the time windows of guidance tasks o and r can be expressed as follows: Figure 10b shows the overlap situation between the time windows of a guidance task with charging demand, and its follow-up tasks.

• Constraints of vehicle sorting
As Fig. 11 is shown, the calculation of assignment inequity index is related to the equity curve and the Lorenz curve.Both of them are formulated on the basis of cumulative distribution of UFMC working time.Therefore, each element of the set of vehicles F = f 1 , f 2 , . . .f n should be sorted in ascending order based on its working time: , ∀l ∈ ∁ F f k

Solution algorithm
In this section, we design corresponding solution algorithm based on characteristics of models mentioned, and construct an overall solution framework for the scheduling problem for UFMCs.As Fig. 12 is shown, the solution process for the scheduling model for UFMCs consists of three steps: preparation of basic data, generation of trajectory on surface, and vehicle resource allocation.
• In the step of basic data preparation, a comprehensive graph model G = (V , E) is generated, where the edges E correspond to various road sections in the scene, including E Apron (the service roads, taxiways and taxi links within the apron), E Taxiway , taxiways in the maneuvering area, as well as E Turn , the various intersections and turns.The vertices V represent the endpoints of the edges, including key points like turning points, entrances/ exits, and crossings.Then, generate a set of candidate routes for each pair of origins and destinations (ODs) on airport surface 26 .• In the step of generation of trajectory on surface, the decision variables R ij , γ ij i ∈ G, j ∈ P are trans- formed into the assigned routes and delivery times in three phases of each guidance task.A function of fitness is formed by integrating the collaborative planning model for surface guidance trajectories and the conflict prediction model based on protection zones.The optimal set of surface guidance trajectories are formed iteratively by the IGWO algorithm.• In the step of vehicle resource allocation, the optimal set serves as the input for the UFMC assignment model.
The model is solved directly using a solver to output the allocation scheme for vehicle resource.
The complexity analysis shows that the scale of feasible solution of the collaborative planning model for surface guidance trajectories is extremely large, with 2 s possible combinations of ONLY the route components.When the number of daily flights exceeds 50, there is likely to be a combinatorial explosion problem.The scale of the decision variable for the UFMC assignment problem is (2ns) , which is a relatively manageable classic ILP problem and does not require a heuristic algorithm additionally.Therefore, to address the collaborative planning model for surface guidance trajectories, an IGWO algorithm based on integer encoding is designed to improve the efficiency of generation for the optimal set of surface trajectories.www.nature.com/scientificreports/

Encoding strategy and initialization
When we apply this algorithm, it is necessary to define the meaning of each dimension in the position vector of a grey wolf individual.Based on the decision variables of the collaborative planning model, targeting guidance tasks for s flights, we construct an integer space with 6s dimensions as shown in Table 2.
Considering that the UFMC scheduling system may avoid conflicts by adjusting routes, shifting time windows, or using a combination of both, conflict-free solution sets would be dispersed at multiple locations, we randomize the wolf pack initially, then integer by rounding to construct the integer space.This serves as the start for iteration and optimization.

Fitness function
The calculation of the fitness function needs to consider the total working time on airport surface, which is based on the selected route, and needs to deduct the situation of potential conflict between the selected route and the variation in delivery time.Therefore, the fitness function f is defined as follows: In this equation, the first term, which relates to the target parking number and selected route, can be obtained by from the pre-computed values of t ij based on the given timetable of flights and information of grey wolves.The second term requires traversing each unit of time for movement on the surface, and ∀node,we need to count the number of times a node falls in different ellipses.

Enhancing global search
In the local search, the slight decimal-level changes in the positions of the wolf pack in the traditional GWO algorithm are not suitable for integer encoding, so we adopt the IGWO algorithm proposed by Yang Z 27 to enhance the prominence and efficiency of global search during the iteration process.The formula for updating the convergence factor is given as follows: Where the improved convergence factor a ′ starts from the default value of 2 and exponentially decays to 2 e with the iteration process.When approximately 69.3 percents of the iteration process is completed, the improved algorithm enters the local search, and the global-to-local search ratio shifts from 1:1 to 7:3.

Improvement for local search
For the local search, we adopt a probabilistic method for updating position 28 : Here, r c is a random number that satisfies U(0, 1) .The formula for updating position based on integer space not only reflects the dynamic process of the wolf pack following the leader, but also ensures the meaningfulness of the operation of updating position for X(t + 1) in the local search. (33) The vector of grey golf X with integer encoding.

IGWO algorithm
In summary, the flow of the IGWO algorithm is shown as following:

Preparation of data
Ezhou Huahu Airport is a specific cargo airport in China, it is also an important experiment field for unmanned ground support vehicles, but there is a lack of overall scheduling for large-scale fleet.Taking Ezhou Huahu Airport for simulation experiment, we investigate the optimization of scheduling for UFMCs during a typical working day in 2030, with 270 flights per day, and assume that runway-to-stand assignments are based on the principle of proximity.Before starting the experiment, we has built a simulation platform.It is divided into air traffic control (ATC) terminal, pilot terminal, and data processing terminal.The interaction between terminals, and the details about hardware and software are shown in Fig. 13.www.nature.com/scientificreports/

Basic routes
Figure 14 shows the schematic diagram of the UFMC's routes.During phases 1 and 3 (operations on service roads), the strategy of separating inner and outer loops is employed.In phase 2 (on-taxiway) operations, the "fixed, one-way, directional, cyclic" principle 29 is followed, where alternative routes are assigned to vehicles, to prevent and alleviate congestion, and to minimize the occurrences of conflict initially.The number of alternative routes, denoted as k, is set to 2. The blue and green arrows in the figure represent Route 1 and Route 2, respectively.
Based on the strategy for route allocation, Ezhou Airport's UFMCs have a set of selectable routes shown in Fig. 15. Figure 15a and b show the departure and arrival scenarios, respectively.In these figures, Route 1 is highlighted in blue (shared with Route 2 in green) during the Guidance phase, while Route 2 is labeled in green.The Dispatch and Recycle phases are marked by the red and cyan lines for Route 1 and 2, respectively.

Parameter settings
According to the actual operation of Ezhou Airport and the requirements of Chinese regulations [20], the parameters of the UFMC are shown in Table 3.The UFMC is allowed to arrive in advance and wait for a maximum of 60s.The time windows for air traffic controllers to adjust the ETD and ETA of flights are 600 s and 300 s, respectively.Therefore, in the case of arrival, the upper limits of the change in delivery times for UFMCs ( γ i1 to γ i3 ) are 60s, 300s, and 60s, respectively; for departure, those are 60s, 60s, and 600s, respectively.

Profiles generation
To illustrate this process, we selected Route 1, along with guarded flights at Stand 111 and 301 as examples, and generated velocity profiles by Algorithm 1. Figure 16 shows the typical velocity profiles.In Fig. 16, the velocity of the GU remains stable at 5 m/s, 8 m/s, and 10 m/s, corresponding to turns, aprons, and straight sections of taxiways in maneuvering areas, respectively.As shown in the second half of each profile, the UFMC adjusts its cruising speed, throttle, and brake, based on a comprehensive modification of the length of each segment and the initial settings, through the application of the Algorithm 1.The velocity is changed when the GU passes through different areas and segments, suggests that the efficiency of the UFMCs scheduling is related to the number of turns, taxiing distance, and timing of crossing the apron area, which confirms the importance of route selection in the trajectory planning module.

Conflicts analysis
An analysis of conflict characteristics before the optimization of scheduling was conducted.Without considering the selection of surface guidance trajectories and adjustments to delivery times, Algorithm 2 predicted 824 instances of conflicts among the trajectories.Figure 17a, b, and c show the characteristics distribution of conflicts spatially, for the west main taxiway, the entire surface, and the east main taxiway, respectively.The characteristic of "high in the north, low in the south" of the distribution of conflicts indicates that, when operating towards the north, the conflicts during arrival are higher than those during departure.On the other hand, due to the coexistence of passenger and cargo parking stands on the west side of Ezhou Airport, the distribution also exhibits the characteristic of "high in the west, low in the east".

Regards to vehicle trajectory planning
After the computation of the vehicle trajectory planning module, the UFMC scheduling system generated a set of conflict-free guidance trajectories, with a total working time of 170,463s.The surface guidance trajectory is mainly planned collaboratively in three aspects: route selection, UFMC's working time, and delivery sequence, achieving the effect of optimization.Among them, the distribution of route selection is shown in Fig. 18.The histogram indicates that Route 1 with shorter travel time remains the top choice for UFMCs.In phases 1 to 3,  Adjustments in UFMC working time are shown in Fig. 19. Figure 19a shows that by extending the working time of UFMCs by 2.85%, the vehicle trajectory planning module achieves the resolution for all the conflicts.Figure 19b and c show that sacrificing operational efficiency is targeted.For arrival guidance tasks with higher conflict frequency, 7.15% of efficiency is sacrificed, while for departure, the change is minimal, with only a 0.66% extension of UFMC working time.
Regarding the result of delivery sequence shown in Table 4, a significant number of guide tasks have been adjusted by fine-tuning the delivery time to avoid conflicts.This has altered the regular delivery sequence,      www.nature.com/scientificreports/

Regards to IGWO Algorithm
To evaluate the convergence speed of the IGWO algorithm, seven different algorithms, including PSO, GWO, WOA 30 , NNA 11,31 , two variations of IGWO algorithm (IGWO-1 and IGWO-2), and the algorithm in this paper (IGWO), were employed in the case study.A population size of 20 and 150 iterations were set, the number of iterations is enough to to make the function variation converge within 1%.The results are showed in Fig. 20.IGWO-1, IGWO-2 and NNA failed to generate any feasible conflict-free solution set for the scale of 270 flights.On the other hand, the remaining algorithms successfully produced feasible solution sets.Among the feasible, GWO exhibited the highest quality but had relatively slower convergence speed.Although PSO and WOA generated the feasible with ideal iterations, their total working times of UFMCs significantly deviated from that of the other algorithms, indicating inefficient operational performance.In contrast, IGWO showed the fastest convergence speed.
Furthermore, in Table 5, quantitative analyses were conducted on four major indicators for handling the scales of 12, 26, 83, 160, and 270 flights per day.In general, IGWO has got the most advantageous indicators (as shown in bold black in the table).WOA is competitive in the convergence speed, but performs poorly in the working time.In contrast, NNA seems to be adapted to small-scale problem.
The performance of GWO, IGWO-1 and IGWO-2 showed improvements compared to the PSO algorithm mentioned in Reference 23 .In the case of 270 flights, the proposed IGWO algorithm achieved an 18.75% improvement in convergence speed, which is consistent with the qualitative analysis.Simultaneously, it incurred a relatively minor loss in operational efficiency, with only 1.76% more total working time compared to GWO.Comparative analysis across different scales revealed that IGWO exhibited strong generalization performance compared to the other algorithms.Meanwhile, IGWO-1 demonstrated potential advantages in solving mediumscale (83 and 160 flight movements) scheduling problems, showing promise for generating more efficient guidance trajectories.www.nature.com/scientificreports/

Regards to vehicle assignment
According to the consultation notice 32 , the estimated number-k of UFMCs in this case study is 9.The results for optimization of UFMC assignment are shown in Fig. 21.The results indicate that each UFMC is assigned 28 to 34 flight tasks, with a standard deviation of 582s in working time.The assignment inequality index is only 0.015731, achieving fair vehicle assignment.Each UFMC only needs to be charged once per day, which is consistent with the total working time of 170,463s mentioned in Section 5.5.When 9 UFMCs are required, the average working time per vehicle is 18,940s, less than half of the safe endurance (q).In fact, for each vehicle, the charging process shown in Fig. 21 needs to be performed only every two days.During the redundant charging time and idle time of UFMCs, they can serve as backups for each other.Considering that Ezhou is a specialized cargo airport, the guidance tasks at night at this airport are relatively active, which is related to the characteristics of distribution of freight flight.

Conclusion
In this study focused on the integrated optimization of scheduling for UFMCs on airport surface, we proposed the concept of GUs and a quantification method for guidance following spacing.We established an integrated UFMCs scheduling model, consisting of trajectory planning and vehicle assignment, considering the demands of precise connecting among guidance tasks and the vehicle charging.Additionally, considering the vehicle dynamics, we designed a high-precision method of guidance trajectory deduction and a conflict prediction model based on protected zones.We also developed an IGWO algorithm fitting the integer encoding.Secondly, we validated the proposed approaches using the case of Huahu Airport in Ezhou.The IGWO demonstrated significant improvements in generating conflict-free feasible solution sets and generalization performance.The IGWO-1 algorithm showed potential in generating high-quality feasible solution sets for medium-scale optimization of scheduling for UFMCs.64.44% and 25% of guidance tasks successfully avoided potential conflicts through fine-tuning timing and detour, sacrificing only 2.85% of efficiency to ensure safety.In the vehicle assignment module, the assignment inequality index was as low as 0.015731, indicating the achievement of collaborative planning and the achievement of balanced allocation of guidance tasks.
Finally, future research will explore tactical optimization strategies for flight delays or unexpected incidents encountered by UFMCs, aiming to enhance the robustness of algorithms and models.

Figure 2 .
Figure 2. The integrated scheduling process for UFMCs.

Figure 3 .
Figure 3. Schematic diagram of a GU.

Figure 5 .
Figure 5. Framework of the scheduling model for UFMCs.

Figure 6 .
Figure 6.Smoothing for segments at turning areas.

Figure 9 .
Figure 9. Framework of the model of UFMC Assignment.

Figure 10 .
Figure 10.Overlap contribution.(a) For time windows of working.(b) For time windows of charging.

( 30 )Figure 11 .
Figure 11.Calculation of assignment inequity index with equity curve and Lorenz curve.

Figure 15 .
Figure 15.Schematic diagram of optional route set for Ezhou Airport.(a) Route set for Departure.(b) Route set for Arrival.

Figure 16 .
Figure 16.The velocity profiles along typical routes.
24.81 %, 21.48 %, and 22.96 % of guidance tasks respectively choose Route 2 to avoid potential conflicts at the cost of detour.

Figure 17 .
Figure 17.Distribution of conflicts before optimization.The southern end of the taxiway corresponds to the right-hand side of the horizontal axis.(a) Conflicts on the west main taxiway.(b) Conflicts on the east main taxiway.(c) Conflicts on airport surface.

Figure 18 .
Figure 18.Histogram for frequency distribution of route selection.

Figure 19 .
Figure 19.Box-plot of UFMC working time.(a) Overall distribution of working time.(b) Distribution of working time for arrival.(c) Distribution of working time for departure.

Figure 21 .
Figure 21.Gantt chart for guidance tasks and charging decisions.The green diamond represents the occurrence of a charging event.

Table 3 .
Parameter settings for the fleet of UFMCs.

Table 4 .
Result of delivery sequence.Significant values are in [bold].

Table 5 .
Quantitative assessment for algorithms.In the same experiment group, the black bold comments mean advantageous indicators.Significant values are in[bold].